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Chiral properties of topological-state loops 
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The angular momentum quantization of chiral gapless modes confined to a circularly shaped 
interface between two different topological phases is investigated. By examining several different 
setups, we show analytically that the angular momentum of the topological modes exhibits a highly 
chiral behavior, and can be coupled to spin and/or valley degrees of freedom, reflecting the nature 
of the interface states. A simple general ID model, valid for arbitrarily shaped loops, is shown 
to predict the corresponding energies and the magnetic moments. These loops can be viewed as 
building blocks for artificial magnets with tunable and highly diverse properties. 
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I. INTRODUCTION 


The appearance of gapless modes at the interface be¬ 
tween insulating phases with distinct topologies is one of 
their quintessential traits, and holds many promises for 
novel device designsj^P— This feature is usually attributed 
to the so-called bulk-edge correspondence. Crudely 
speaking, in order to “unwind” one band structure into 
another one must close the gap at some point, render¬ 
ing protected localized states, that usually display chi¬ 
ral propertiesi^ The prime example of this is the Quan¬ 
tum spin Hall (QSH) effect. There, counter-propagating 
modes belonging to opposite spins will emerge within 
the bulk band gap at the boundary between two time- 
reversal symmetric insulators with different Z 2 topologi¬ 
cal invariants i^^— 

While the manifestation of chirality in electrical trans¬ 
port has been thoroughly researched,^— in this paper 
we will examine in detail what happens when the chiral 
states are forced into a loop delineating two topologically 
distinct insulating domains. This is somewhat different 
and more general than the case studied in Ref. [l^, where 
the impact of circular vacuum domains in Z 2 topological 
insulators on transport was studied. In particular we will 
show that, depending on the particular topology and the 
type of material on either side of the loop, the chirality 
manifests as a spin and/or valley coupling between the 
states bound to the orbit and their angular momentum, 
making these loops potential building blocks for hybrid 
magnetic structures. More importantly, we will also show 
that a general ID model captures all of the physics in¬ 
volved, even for noncircular loops, and discuss the limited 
impact of intervalley scattering. Our work is in partic¬ 
ular focused on group IV honeycomb monolayers, which 
we analyze within a continuum approximation by em¬ 
ploying the Dirac equation and the 7r-band tight-binding 
(TB) model. 

The paper is organized as follows. In section [TT] the 
analytical solutions for circular interfaces are studied us¬ 
ing the Dirac equation. In section m the tight-binding 
method is used to examine the impact of the magnetic 
fields on the topological states for varying loop geome¬ 


tries, as well as to establish the validity of a general 
formula for energy levels and magnetic moments of the 
states bound to the loops. Finally, in section HVl we dis¬ 
cuss the practical consequences of our proposal and sum¬ 
marize our results. 


II. THE CONTINUUM APPROACH: 
CIRCULAR INTERFACES 


We start by introducing our general model, which is 
the Dirac equation valid for honeycomb-lattice materials 


H = hvF {rkj^ax + kyay) + STAKMcrz+TAHcr2 + Aspaz. 

( 1 ) 

Here r = -|-I (r = —1) denotes the K {K') valley and s = 
-l-I (s = —1) labels spin up (spin down). Akm is Kane 
and Mele’s phenomenological term modeling spin-orbit 
coupling,— while Ah was introduced by Haldane as a toy 
model for a non-trivial band structure in the absence of 
magnetic fieldfiS inducing the quantum anomalous Hall 
(QAH) phase. Finally, Asp is the staggered potential 
breaking the inversion symmetry, and opening a trivial 
band gap unlike the previous two. 

An interface between any of the aforementioned insu¬ 
lating phases will support boundary modes. This can be 
shown by solving Eq. ([1]) for the geometry depicted in 
the upper row of Fig. [TJ Note that in the domains / 
and II only a single mass term of the form Aij^az is 
used in Eq. O- Eor instance, in case (c) A/ = stAkm 
and Ajj = Asp. Then the evanescent modes, decaying 
exponentially away from the interface read 
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( 2 ) 


and we are interested in solutions E = Ahvpkx with 
i^yi/ii = \^i/ii\ /^tip. Matching the wavefunctions at 
the boundary leads to the dispersion relations shown in 
the middle row of Eig. [TJ These demonstrate chirality of 
the topological modes bound to the interfacej^ Unlike 
the cases in Eigs. [T](c) and (d), the first two cases are 
spin degenerate. 





2 



FIG. 1. (a) - (d) The considered set of interfaces between in¬ 
sulating phases. The upper row shows the parameters related 
to each region, while the middle row depicts the corresponding 
dispersion of the emerging modes. The bottom row illustrates 
the hypothesis for the persistent currents in the closed loops 
inferred from the linear dispersion relations. Red (blue) lines 
depict the states belonging to the K {K') valley, while the 
spin is denoted by the projection of the arrows. 


TABLE I. Summary of the numerical solutions. The symbol 
“x” means that there are no allowed solutions, while the al¬ 
lowed modes are characterized by the corresponding angular 
momentum numbers j. 



t K 

t N' 

IK 

IK' 

(a) 

j>0 

i < 0 

j > 0 

j <0 

(b) 

X 

i < 0 

X 

j <0 

(c) 

X 

j <0 

i > 0 

X 

(d) 

X 

X 

j <0 

j < 0 


We now present the main premise of our paper. Instead 
of straight interfaces, we consider circular ones, depicted 
in the bottom row of Fig. [T] The chirality of the bound¬ 
ary modes suggests that the angular momentum of the 
states in these loops must obey certain selection criteria 
depending on the particular phases in regions I and II. 
For instance, in the cases (b) and (d) one would expect 
to see that only negative angular momentum valley and 
spin filtered states are allowed, respectively 

To explore this explicitly, we solve the problem in polar 
coordinates, in which the Hamiltonian ([ij reads 


A//// TT 

7r+ -A//// 


( 3 ) 


where {—idr ± ^d^). Because of the 

radial symmetry the total angular momentum Jz = Lz + 
is a good quantum number denoted by and 

therefore the solution can be sought in the form 




XA{r) 


( 4 ) 


The coupled system of equations reduces to the differen¬ 
tial equation 


dr'^ 


dr 


(5) 

Having in mind that we are interested in states lying 
within the bulk band gap > Q, so that the 

solutions are the modified Bessel functions Ij-r/ 2 {'i’) and 
Ej-r/ 2 (j')^ which are exponentially divergent at r —> oo 
and r = 0, respectively. The radially bound modes read 


and 


-f j-r /2 {z ijr)) 
rd>V^f^r 

^ i 


Ai+E 


, ( 6 ) 


'i’ll = e 




Kj-r /2 (ziijr)) 


ire 


r^VMEl! 


An+E 


Kj+rl 2 {zii{r)) 


where zj/u = 




( 7 ) 


V/// 


-E^ 


interface r = R results in the eigenvalue equation 


r. Matching vk/ and i’ll at the 


/A^ — 

Au + E 

/^2 _ j ^2 

+ a/+f; 

which we solve numerically. 

Our numerical results are summarized in Table □ for 
four cases shown in Fig. [TJ For those spin and valley fla¬ 
vors for which there exists a solution for a given interface 
we give the allowed angular momentum values, while the 
cases for which no solution is found are indicated by “ x ”. 
Only results for positive energies are given, since the neg¬ 
ative energy solutions are the same, albeit with opposite 
values of j. By inspecting Table H] and comparing it with 
Fig. [1] one can see that the numerical results agree with 
the expectations deduced from the chirality of the topo¬ 
logical modes. Namely, the angular momentum of the 
states bound to the loops not only display the expected 
orientation but are also: (a) valley-coupled, (b) valley- 
filtered, (c) spin-valley-filtered, and (d) spin-filtered. 

In Fig. [2]the solid black lines are the numerically com¬ 
puted energies of the bound states for the radial interface 
corresponding to Fig. [T](a) for spin up in the K valley. 
Note that for all the modes that are allowed, the set of 
energies is the same (only the angular momentum num¬ 
ber j varies), such that it suffices to show the spectrum 
for only one particular case. One can see that the larger 
loops can support states with larger j. Note also that, 
as a rule, there exists only one state for a given j. This 
is due to the exponential localization at the loop, which 
prevents oscillatory behavior in the radial direction and 
the emergence of a corresponding degree of freedom. 

Interestingly, the energy curves follow a 1/R behav¬ 
ior with very high precision, a property also seen in 
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R (nm) 

FIG. 2. Variation of the numerically computed energy lev¬ 
els with radius R (solid black lines), and the results from 
semiclassical quantization (dashed red lines) for the loop con¬ 
sidered in Fig. [T](a). Here Asp = 0.5 eV. 


states bound to the holes punched through topological 
insulatorsi^ It is relatively easy to obtain a good ana¬ 
lytical fit of the energy levels using the following reason¬ 
ing. Given the exponential localization of the modes near 
r = R, the loops effectively become ID closed lines carry¬ 
ing states with pseudospin degree of freedom, described 
by ■)/). Then, by requiring single-valuedness of "0 upon 
one full revolution we must have ipi'P + 27r) = 
where the contour is oriented counterclockwise. The mi¬ 
nus sign comes from the spinorial nature of the state, and 
it forces the phase factor accumulated along the trajec¬ 
tory 6 = sgn{j)pL/h {L = 2 'kR is the loop length) to be 
quantized in half integer multiples of Denoting this 
half integer number by j, and recalling the dispersion 
of the topological modes E = vpp, yields the following 
expression 


E=\j\ 


hvp 

~R' 


(9) 


for the energy levels. They are shown in Fig. [2]by dashed 
red lines, and are in excellent agreement with the numer¬ 
ical results. Only small deviations are found for loops 
with small radius. This is to be expected, given that for 
small loops overlap of the wavefunction must occur, per¬ 
mitting tunneling of the modes between sections of the 
loop. 


III. THE TIGHT-BINDING METHOD: 
GENERALIZATION FOR ARBITRARILY 
SHAPED LOOPS 


Having analytically elucidated the chiral nature of the 
angular momenta, we now employ the TB method to 
study the behavior of the loops in a magnetic field. On 


the one hand, this is important because it can offer in¬ 
sights into intervalley scattering. Indeed, the most crit¬ 
ical case is the setup (a), since valley flipping can cause 
back-reflection and mixing of the states, something not 
captured within the continuum approach. All the other 
loops are more robust, since either there are no counter- 
propagating states ((b) and (d)), or backscattering re¬ 
quires the simultaneous flipping of the valley and spin 
indices (setup (c)). On the other hand, the TB method 
will help us to resolve important details of the magnetic 
moments originating from the loops, even for noncircu¬ 
lar geometries. The TB Hamiltonian that we employed 
reads 

H — -t 

±± — 1/ y ^ o l 0^0777,. 

{n,m) {{n,m)) 

( 10 ) 

We use t = 2.7 eV for the hopping between nearest neigh¬ 
bor orbitals, while the second term is the SOC. Note 
that Vnm = +1 {vnm = —1) if an electron makes a right 
(left) turn at the intermediate atom when hopping from 
site m to site n. The Peierls term (pnm = ~'E fv" 
accounts for the phase the electron acquires while travel¬ 
ing in the presence of the magnetic field. The staggered 
potential is simply added along the main diagonal of the 
Hamiltonian, like in the continuum model. 

Our calculations were performed for a loop of radius 
i? = 10 nm contained within a larger circular quantum 
dot, with a 15 nm radius, so that a proper decay of the 
topological modes is ensured. The results are shown by 
solid lines in Figs. Oja) and (b) for the loop setups de¬ 
picted in Figs. (Ha) and (c), respectively, as a function of 
the magnetic flux through the loop <i> = BR^tt in units of 
the flux quantum 3>o = hje. For setup (a) the results are 
shown in black, since the TB model cannot resolve the 
two valleys, while for setup (c) (Fig. [3Kb)) the valleys can 
be distinguished due to the additional spin polarization. 
In Fig. ISKa), for i? = 0 the majority of the levels are de¬ 
generate indicating the preservation of the valley index. 
Only a handful of levels (for instance those around 0.1 
eV) are offset by an equally small amount in the opposite 
directions with respect to the continuum solutions. This 
behavior indicates that such states are susceptible to in¬ 
tervalley scattering, causing bonding and anti-bonding 
states. As the magnetic field is turned on, the remain¬ 
ing valley degeneracy is lifted due to the opposite mag¬ 
netic moments (see below). Moreover, an anticrossing 
behavior occurs at energies where intervalley scattering 
is prominent, resulting in an Aharonov-Bohm-like pat¬ 
tern of the curves. On the other hand, for the setup (c), 
such scattering is impossible in the TB picture due to 
the absence of spin-flipping terms in the Hamiltonian, 
as demonstrated by Fig. ISKb). Consequently, the spin 
up in the K' valley (solid blue curve), with a magnetic 
moment oriented along the magnetic field, and the spin 
down in the K valley (solid red curve) with the opposite 
orientation of the magnetic moment are independent of 
each other, and experience a rigid, linear, decrease and 
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FIG. 3. The magnetic field dependence of the energy spec¬ 
trum for the loops depicted in (a) Fig. [2a) and (b) Fig. [2c). 
We used Asp = Akm ~ 0.5 eV. Solid curves depict the TB 
results, while the dashed lines are given by Eq. dill) . Red 
(blue) depicts the states belonging to the K {K') valley. 


an increase in energy, respectively. 

In order to capture the behavior of the topological 
modes in the presence of a magnetic field, we again con¬ 
sider the phase change upon a single rotation around the 
loop. Now, due to the magnetic field an additional phase 
factor of the form (assuming B is constant inside the 
loop) — I ^ A • dl = —j-BS appears, where S is the loop 
area. As before, 6 must be a half integer multiple of 27r, 
yielding 


E = |j| ^ + sgn{j)^SB. 


( 11 ) 


These levels are depicted by dashed curves in Fig. [3] and 
show excellent agreement with the TB results, with small 
deviations at very large magnetic fields.— Note that the 
change in energy due to magnetic field can be captured 
entirely by the classical formula for the energy of a closed 
current contour given by —fiB. Besides the area of the 
loop, the magnetic moment of a classical contour fj, = IS 
depends on the current due to one topological mode I = 
sgn(j)(—e)/T, where T = sgn{j){—e)vF/L is the period, 
and matches the one in expression Eq. (22 exactly. 

Most importantly, note that Eq. 22 is written in 
a general form, without specifying the precise shape of 
the loop, since it depends only on its length and area. 
In fact, we have verified numerically by using the TB 
method that the above formula accurately predicts the 
energies and magnetic moments of the states in loops 
with arbitrary shape. In particular, in Fig [2 we show the 


FIG. 4. The same as in Fig. [3l but for an elliptic loop ((a) 
and (b)), hexagonal zigzag loop ((c) and (d)), and hexagonal 
armchair loop ((e) and (f)). 


results for the same two setups as in Fig. [3l but for a 
elliptic, hexagonal zigzag and hexagonal armchair loops. 
Note that j does not denote the total angular momentum 
now, due to the lack of rotational symmetry. However, it 
is still required to be a half-integer, while its sign reflects 
the orientation of the chiral current. For the elliptic loop 
(Figs. 12a) and (b)) we have used i?i = 5 nm, i ?2 = 15 
nm, S = i?ii?27r, and for the loop length we have used 
the approximation of the form 

L = 7r(i?i-hi?2) 1 1 +-(12) 

where d = (i?i — i? 2 )^ / (.Ri + The hexagonal 

loops are distinguished by the number of zigzag {Nzg) 
or armchair (Nac) segments along one edge. The loop 
length {L = 61) and area {S = 6P-\/3/4), depend on 
the length of one side of the hexagon Izg = NzgV^g 
{IAG = {^Nag — 2)a), where a is the nearest neighbor 
distance. The results depicted in Figs. [2c) and (d), and 
(e) and (f) are for Nzg = 42 and Nag = 25, respectively. 
All the loops are surrounded by a region large enough to 
ensure the decay of the topological modes. 

On the one hand, for setup (c) (where counterpropa- 
gating states belong to opposite spins, lower panels), it 
can be seen that regardless of the shape, the persistent 
current is robust, and that Eq. 22 predicts the energy 
levels with great accuracy. The very small mismatch that 
appears we attribute to the fact that the sharp turns and 
corners can effectively reduce the loop length and area 
due to tunneling, and hence cause some deviation from 
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the expected behavior. On the other hand, the spectrum 
of setup (a) (upper panels) shows that intervalley scatter¬ 
ing is geometry dependent, as one would expect, although 
Eq. (HB) still provides a very good approximation for the 
results in general. While elliptic loops display stronger 
intervalley scattering and subsequent hybridization than 
the circular loops, inter valley scattering in zigzag hexag¬ 
onal loops is practically nonexistent. Finally, in arm¬ 
chair hexagonal loops a gap is opened for low-energy chi¬ 
ral states, a situation resembling the appearance of mass 
barriers in graphene quantum rings Remarkably, away 
from this energy region a behavior indicating chirality 
preservation is observed, although the valley degree of 
freedom is ill-defined in this structure. Therefore, to a 
good approximation one can say that the electronic and 
magnetic properties of these loops are largely indepen¬ 
dent on the precise details of the edges near the interface. 
This is unlike the case in regular graphene nanorings, 
where geometry and boundary of the rings can hinder 
the fidelity of the simplified ID models;^ with the main 
reason for this contrast being the exponential localization 
of the topological states. 

IV. DISCUSSION AND THE CONCLUSION 

Now we consider the possible applications of these ef¬ 
fects. By controlling the valley population separately in 
the setup (a) for example (in loops with well preserved 
valley index), one could in turn control the total mag¬ 
netization of the loop, since all orbits within a given 
valley have the same orientation of the persistent cur¬ 
rent. In setup (b) the K valley does not host topolog¬ 
ical boundary modes, unlike the K' valley, where the 
persistent current has the same orientation regardless of 
spin. Hence, these particular loops should be inherently 
magnetic. Alternatively, in setup (c) one can control 
the total magnetization by adjusting the spin or spin- 
valley populations)^ Finally the loops in setup (d) not 
only display nonzero magnetism, as those in setup (b), 
they have the additional benefit that the allowed topolog¬ 
ical modes have their spins aligned, which can contribute 
to the overall magnetic moment. Therefore, these loops 
have a tremendous potential for the design of novel mag¬ 
netic structures and devices. One can imagine an array 
of these loops laid out in a pattern designed to suit a 
particular application. Moreover, since the states racing 
along these loops are exponentially localized in the radial 
direction the array can be densely packed, enabling the 
creation of hybrid magnets. 

Before we conclude, let us reflect on the practical is¬ 
sues concerning the real materials that are required to 
create these structures. On the one hand, it is well 
known that inversion symmetry can be broken by lay¬ 
ering graphene with boron-nitride, which has a natu¬ 
rally occurring staggered potential since it’s sublattices 
are composed out of different atomic speciesi^ Moreover, 
in-plane heterostructures of graphene and boron-nitride 


have been demonstrated in recent experiments In 
particular, in Ref. the authors presented a matrix ar¬ 
ray of circular graphene dots embedded in boron-nitride. 
Since the resulting structure was transferable, one can 
imagine layering it with a substrate which would enhance 
SOC in graphene, such as WS 2 for instance^ thereby 
creating a macroscopic pattern of chiral closed-contour 
currents of the type shown in Fig. me). 

On the other hand, given the buckled lattice structure 
of silicene, germanene and stanene^^— breaking of the 
inversion symmetry can be realized by external gates, 
while relatively heavy constituent atoms ensure nonzero 
SOC. Patterning gates on top of them could likewise in¬ 
duce chiral persistent currents described in this paper, 
making them detectable for instance by magnetic force 
microscopy. Additionally, there are theoretical propos¬ 
als for inducing the QAH effect in graphene and silicene, 
although experimental validation is still lackingi^^— 

As we have shown, the larger the bulk band gap and 
the loop radius, the larger the magnetization will be. 
However, large loops and inherent spin and valley re¬ 
laxation mechanisms will inevitably lead to scattering, 
spoiling the effect to a certain degree. Nevertheless, no 
amount of scattering can change the orientation of the 
persistent current, which is related to the topology of 
the surrounding band structures, and the emptied states 
should quickly become repopulated. Note that, while 
the interplay of orbital magnetism emerging in insulat¬ 
ing Dirac systems)^^— and the magnetism of topological- 
state loops remains to be explored, the main obstacle for 
realizing this effect is most likely related to the size of the 
bulk band gaps. This can not only restrict the amount 
of persistent current supported by the loops, but more 
importantly limit the temperature range of applicability. 

To summarize, we investigated the behavior of elec¬ 
tronic states bound to looped interfaces between insulat¬ 
ing phases of distinct topologies. These structures can 
be a source of magnetic moments due to persistent cur¬ 
rents flowing along the interfaces. The chirality of the 
states bound to these loops manifests as a chirality in 
the angular momentum and the corresponding magnetic 
moments. We have shown that a simple ID model cap¬ 
tures both the qualitative and the quantitative behavior 
of the electrons. By calculating energy levels and mag¬ 
netic moments in elliptic, hexagonal zigzag, and hexago¬ 
nal armchair loops, it was shown that the model’s validity 
extends to noncircular loops as well. Because of the inti¬ 
mate link between the magnetic moments and spin and 
valley degrees of freedom, novel magnetic structures and 
devices can be envisioned. 
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